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Stationary axisymmetric perfect fluid space-times are investigated using the curvature 
description of geometries. Attention is focused on space-times with a vanishing electric 
part of the Weyl tensor. It is shown that the only incompressible axistationary magnetic 
perfect fluid is the interior Schwarzschild solution. The existence of a rigidly rotating 
perfect fluid, generalizing the interior Schwarzschild metric is proven. Theorems are 
stated on Petrov types and electric/magnetic Weyl tensors. 

1 Introduction 

A well-known property of the static, incompressible interior Schwarzschild space-time is con- 
formal flatness, i.e., the Weyl tensor vanishing: Cijki — 0. The Weyl tensor can be split 
into two parts, the electric and magnetic parts. The first part represents the tidal forces of 
the curvature, while the second part (which has no Newtonian analogue) is generated by the 
vorticity, shear, etc. of the congruence to which the electric and magnetic parts are related. 
Therefore, rotating matter can be viewed as partly responsible for the existence of the mag- 
netic part of the curvature Q. A simple example of a space-time containing rotating matter 
is expected to be a perfect fluid with a purely magnetic Weyl tensor. 

Little is known about space-times with a magnetic curvature. Arianrhod et al. ^ have 
obtained an example of these space-times, albeit with an unknown equation of state. Another 
solution, suggested in is, however, shown to be not purely magnetic by Mcintosh 

To get some more insight on the subject, in this paper we investigate perfect fluid space- 
times with axistationary symmetry. We confine our attention to circularly rotating matter, 
i.e., the four- velocity lying in the plane of the two Killing vectors. Thus we may choose the 
line element in its Papapetrou form, 

ds^ = f{At + LO A^f - /-I [e^^ (dr2 + dz^) + d^^] . (1) 

We choose the canonical frame 

0" = (dt + tjd(p) 6/i = /-i/Vdr 



02 = /-i/2eMz 03^/-i/2pd(^ 



(2) 



and project all quantities onto this. 

A space-time is locally completely determined by the set 



S = {Rijki , l%k ,eix", nij } , (3) 
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where Rijke are the coraponents of the Riemann tensor, are the rotation coefficients, 
Bix" the gradients of the essential coordinates, and rjij is the constant metric in the chosen 
frame^ ^. The set S may be calculated using the frame (|^) and has a form characteristic 
of the metric (|^). We still have some freedom in performing Lorentz transformations of the 
frame (|^), e.g., rotations in the (l,2)-plane (see Sec. ^). These transformations will change 
the appearance of the components of, say, the Riemann tensor, but will not introduce new 
non-zero components of this tensor. This freedom can be exploited to make the set S as 
simple as possible. There is also a boost freedom, which does not introduce any new non-zero 
components of the tensors in the set S. This can be used to put the frame (||) in a comoving 
form, thus giving a natural interpretation of (some of) the Ricci rotation coefficients in terms 
of fluid variables (such as vorticity). 

If we start with a set S, instead of generating the set S from a metric satisfying Einstein's 
field equations, we need to impose certain restrictions on the elements of the set in order to 
guarantee that S can be generated by a line-element that satisfy Cartan's equations under 
the restriction of Einstein's field equations. In general, these conditions are part of the 
commutator equations, Ricci equations, and Bianchi equations ^. 

The approach taken in this paper proceeds in the following steps: 

1. Investigate the form of S in the Lorentz frame (||) for the metric (|l|). 

2. Analyze the properties of this S under Lorentz transformations. 

3. Using the knowledge obtained in the steps above, make an ansatz for the set S. 

4. Subject the elements of this Ansatz to the necessary and sufficient equations (in general, 
parts of the commutator equations, Ricci equations, and Bianchi equations). 

5. Use the fact that part of the variables in S are the components of the Weyl tensor to 
investigate magnetic space-times. 

In Sec. H we derive the equations describing an axisymmetric and stationary perfect fluid 
space-time. These are the commutators acting on the "essential coordinates" (on which all 
quantities in S depend) and the Ricci equations. The Bianchi equations are not needed, 
but in practice they are very useful when there are restrictions on the Weyl tensor. All 
these equations are first order partial differential equations in terms of two variables for the 
quantities in the set S. Also, S is chosen as to be compatible with the metric (|^) in a Lorentz 
transformed version of the frame (^ . 

In Sec. ^ we investigate the restrictions of Sec. ^ in the case of a purely magnetic Weyl 
tensor. Specifically, we examine what space-times are allowed under the further restriction of 
incompressible fiuid (constant energy density) or shear- free motion of the fluid (rigid rotation). 

Section ^ describes how to obtain the full metric of the space-time. We focus on the 
physically important shear free space-times. We show that the line-element can be written in 
terms of quadratures once the elements of the set S are determined by the equations set up 
in Sec. |. 

In Sec. |, we consider various Petrov types of circularly rotating perfect fluids. We prove 
that type N fluids have the unphysical equation of state /i + p = 0, and that rigidly rotating 
purely electric fluids are of Petrov type D. 

^Note that the analysis would be possible in any constant frame, not just the Lorentz frame. We choose 
this particular frame since it is related, in a simple fashion, to the physical observables of the fluid. 
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2 Equations 



Our space-time possesses two Killing vectors, one time- like and one space- like, which commute. 
We assume that the coordinates are adapted to the Killing vectors, such that all quantities 
in the set S are dependent on two coordinates {x"} = {x,y}. 

The source is a perfect fluid with the energy-momentum tensor 



(4) 



We choose a frame = {eiX^){d/dx^) with dual a;* = lj'^ da;^, i ~ 0, 3, such that 

= - (,^1)2 _ (,^2)2 _ (^3)2 

Co = M , 



(5) 
(6) 



i.e., the frame is comoving and (.9^) = {flij) = diag(l, —1,-1,-1). In terms of the frame ^ 
related to the metric (|l|), the choice of a comoving frame means that we have carried out a 
boost in the (0,3)-plane: 



(7) 
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or explicitly 



6° 



9" = coshx/^^^(di-f wd(^) +sinhx/"^^^(?d(p 
9^ = coahxf g dip + sinhxf^^'^ [dt + uj dip) 



(8) 
(9) 



where x = xi^iU)- In. the case of a rigidly rotating fluid, the boost parameter x may always 
be put to zero. With this, the contravariant vectors Bq and 63 become 



eo 



63 



1 r 

Q 

1 r 



coshx/"^/^£') dt ~ sinhx/^/^9>. 



fsinhx/^/'^ 



fcoshx/^/^w sinhx/^^/^fi*) dt -f coshx/^/^5, 



(10) 
(11) 



We define the non-zero components of the Riemann tensor of the circularly rotating fluid 
in the frame given by Eqs. (|^) and (^ as 



^0101 

^0113 
^0202 
^0223 
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^1313 
-R2323 



Si-i(Ai + 3p) i?o 



^0123 — 



■E^2-^(M + 3p) i?0213 



^3 

—H\ — H2 
-E2 

-El - ^p , 
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E3 

Hi 

-H2 

-Ei-E2-\[p + ip) 

E1+E2- i/i 

E3 



where Ea and Ha, A = 1,2,3, are components of the electric part Eij = C, 



ikjel 



,k„,l 



and 



magnetic part Hij = *Cikjiu^u^ of the Weyl tensor respectively. The rotation coefficients 
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7jfe = ijj\.t,{ejxi-'){ekx'') are 
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where is the acceleration, ai the shear, and Wi the vorticity of the fluid. 
We stih have the freedom to perform rotations 
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(12) 



in the (l,2)-plane (here = 'd{x,y))^ This rotation wiU not introduce any new non-zero 
components in the Riemann tensor or the rotation coefficients (as any other rotation would 
do), although the components of the respective geometric objects will mix. 

Parts of the Ricci equations and Bianchi equations are solved by the choice Bqx = Cqj/ = 
e^x = e^y = 0. Using the reduced number of coordinate gradients, we may write 



ei = £.idx + vidy 
62 = C2<9x + y2dy , 

where = eix, Vi = eiy. 

The commutator equations are generated by 

We will refer to a specific commutator by the index pair (i, j) 
tator on x and y, we obtain the 



commutator equations 



ei6 - 626 = -aiCi + "2C2 
eiV2 - 62^1 = -aiVi + a2V2 



(13) 
(14) 



(15) 

Acting with the (1, 2) commu- 



(16) 
(17) 



Writing Bq and 63 as linear combinations of dt and 9^, from the (0,3) commutators acting 
on x and y follows that 



71 
72 



(Tl + UJl 
(J2+UJ2 , 



(18) 
(19) 



respectively. All the remaining commutator equations are identically satisfied when acting on 
the essential coordinates x and y. 

^ During this transformation, the rotation coefficients transform as 

7> = A\Ar^k"Ym„ + (e„A',)A/Afc" , 
where A'^ is the transformation matrix and A*f A/ = 5'.. 
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From the relation R%ij = eij%j — ejj%i + 27fcb[i7'''^j] + 2j%k-f\ij] we obtain the 
Ricci equations 

Si - + 3p) + al + 02^1 - (0-2 + t^2)(3cr2 - UJ2) (20) 

E2 - + 3p) + al + aia2 - (cri + ti;i)(3a-i - uJi) (21) 

-B3 + aia2 - a2Q!2 + (3cri - c^i)(o-2 + ^2) (22) 

E3 + 0201 — aiai + (3(T2 — tt^2)(o'i + ivi) (23) 

-iJ2 - 2a2W2 + ai(o'2 + t^2) - /3i(o'i - ti^i) + /32(cr2 - UJ2) (24) 

-iJi - 2aiWi + q;2(o'i + t^i) - /32(o'2 - ^^2) + A(cri - ^1) (25) 

H3 + 2a2UJi - a2(o'2 + UJ2) - 2/32cri (26) 

H3 + 2aiLU2 - ai{ai + uji) - 2l3ia2 (27) 

-E2 - + ai/32 - /?! - (ct2 - 3tJ2)((T2 + UJ2) (28) 

-^1 - i^i + a2/3i - I3l - (ai - 3iUi){ai + Wi) (29) 

i?3 - a2/32 - /?i/?2 + (di - 3ci;i)(cr2 + W2) (30) 

E3 - ai(3i - /32/3i + (0-2 - 3w2)(o'i + wi) (31) 

and 

ei((Ti - uji) + 62(0-2 - UJ2) = -Hi - H2 + ai(o-2 - ^2) + "2(0-1 - wi) (32) 

6201 + eiQ;2 = -i?i - i?2 + + Q^i + ^2 (33) 

£^1 + £^2 = -^(a* + 3p) - aiPi - 02/32 + (^i + uif + (0-2 + <jJ2f (34) 

Hi+H2 = -(ai +/3i)(cti +cji) - (02 + /32)(o-2 + ^^2) ■ (35) 

The Bianchi equations Rij[kt,m] — necessary equations (since the space-time is 

■without isotropy but they may be used as "help" equations (although all information 

in them are contained in the lo-wer-derivative stage of the Ricci equations for stationary 
axisymmetric space-times). Thus we have the 
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twice contracted Bianchi equations T^-'-j 



eip = ai{p+p) (36) 
62P ^ a2{pL+p) , (37) 



and the 

Bianchi equations 



61 (£1 + E2) = (ai - 2/3i)£i + (2ai ~ Pi)E2 - (02 + /32)-B3 

+ (ai + LUi){Hi - H2) + 2(ct2 + W2)i?3 - ^eiAi (38) 
62(^2 +Ei) = {a2- 2/32 )£2 + (2a2 - P2)Ei - (ai + l3i)E3 

+ (02 + LJ2){H2 - Hi) + 2{gi + L0i)H3 - ie2M (39) 
61^1 + 62^3 = (02 - 2(3i)Ei - {a2 + Pi)E2 + (2ai - /32)£3 

+ (cri + 3tJi)iJi + 2o-iiJ2 - (0-2 - 3cj2)-ff3 - (40) 
62^2 + ei£;3 = (ai - 2P2)E2 - [ai + (32)Ei + (2^2 - /3i)-B3 

+ ((72 + 3w2)i?2 + 2o-2i^i - (0-1 - 3wi)ff3 - |e2M (41) 
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and 



ei(iJi + H2) = (ai - 2(3{)H^ + (2ai - /?i)F2 - (02 + l32)H:i 

+ (o-i +cji)(£;2 - El) - 2{a2+LU2)E3 + +tJi)(Af +p) (42) 
62(^2 + i/i) = (02 - 2/32)i/2 + (2a2 - /32)ifi - (ai + /3i)i?3 

+ ((72 +c^2)(^i - E2) - 2(cri +cji)£;3 + i(cr2 +t^2)(Ai+P) (43) 
eii/i + e2i/3 = («2 - 2/3i)i/i - (aa + /3i)i/2 + (2ai - ^1)11^ 



Once more, we stress the fact that in order to find a solution to Einstein's equations, we 
need to solve Eqs. (|T|), (|T|), and (|ol)-(|5l). The use of the Bianchi equations (|3|)-(||) is 
optional when no further symmetry restrictions are imposed. 



3 Magnetic space-times 

For space-times with a purely magnetic Weyl tensor, we set Ea = 0. 



3. 1 Incompressible fluids 

Incompressible fluids are defined by ei^ = 62^ = 0. We use the rotation (|l^ ) to set = 0. 
We may then prove the following: 

Theorem 1. The only incompressible axistationary perfect fluid with a pure magnetic Weyl 
tensor is the interior Schwarzschild solution. 

Proof. From Eqs. (^8|) (^l|) we get four algebraic equations from which Hi — H2 and 
ai + u!i — a2 + 1^2 = 0. But Eq. (js^) then results in Hi — H2 — 0, i.e., the space-time 
is conformally flat. According to a theorem by CoUinson 0, the only conformally flat axista- 
tionary perfect fluid space-time is the interior Schwarzschild solution. □ 



3. 2 Rigidly rotating fluids 

When the fluid is rigid, we have (Ti = (T2 = 0. We set o^i = by using the transformation 
(p^. Then the Bianchi equations (p8|)-(|4l|) give the constraints 



62^2 + eiH^ 



(cTi -I- iuuijEi - 2(71 £^2 + (0-2 - 3w2)i?3 + +p) 

(ai - 2/32)i?2 - (ai + /32)i?i + (2a2 - /3i)i?3 

((T2 + 'iL02)E2 — 2a2Ei + (di — 3uJi)E3 + u;2{^^ + p) ■ 



(44) 



(45) 



0J2H3 = 0, Lj2{2Hi + H2) = . 



(46) 



If UJ2 = 0, the fluid is static and Eqs. ( p8| ) and ( p9| ) imply ei/i = 62/1 = 0, while Eqs. (p^, 
(p6|), and ( p5| ) gives Hi = H2 = H3 = 0. Thus, this is the interior Schwarzschild solution. 
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If i/3 = and H2 = 



—2Hi, we obtain a consistent subsystem of equations 



eip = Bill = ei/32 = 6102 = (47) 

eafla =03- ^(M + 3p) (48) 

e2|32 = -\^^- Pi (49) 

e2M==-3(a2 + /32)(yU + 3p + 602/32) (50) 

e2P = a2(/i + p) (51) 

ei/3i = - /32 (52) 

e2/5i = -/5i/92 , (53) 

where 6^ = /3| — 3a2/32^^M^|p- (For the interior Schwarzschild solution 6^ > 0.) This system 
is characterized by wi — ai = a2 = 0, cti — —P2, Hi = {02+^2)^21 anda;| = (/x+3p)/6+a2/32. 
We choose coordinates {x, y} such that ei — Q{d/dx), where we assume that C, = C{y) depends 
only on y, and 62 = ^{d/dx) +v{d/dy). With this choice of coordinates, all quantities except 
Pi become independent of x. We still have the freedom to make the coordinate transformations 

X = hiiy)x + hoiy) (54) 

y = y{y) , (55) 

where hi{y) and ho(v) are some x-independent functions. It follows from the (1,2) commu- 
tator on yicf. Eq. (|l7|)), that v depends only on y. Acting with the commutator (1,2) on x 
{cf. Eq. (©), we find that 

^T- = CP2 ■ (56) 



' dx dy 



Integrating, we obtain 



C=(/32 + ^|)x + eo, (57) 

where ^0 is a function depending only on the coordinate y. Performing a transformation (^^, 
^ becomes 

^^^hi{y)+v[x + j ■ (58) 

Substituting ([s^), we can achieve f = by suitably choosing hi{y) and hoii/). This means 
that we can use a coordinate system in which 

-i=C^. ^2 = ^.^, (59) 

and both C and v depend only on y. Since b does not depend on x either, we can integrate 
53) to obtain 



Pi — —5 tan 



(60) 



where xq is some function of y. Substituting ^ and Pi in Eq. ( |53| ) and using 626 — —P2b we get 
that Xq must be a constant. Since we only used the derivative of ho{y) in the transformation 
we still may add a constant to the coordinate x, and we use this freedom to set xq — 0. 
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The remaining equations are 



v^ = al^l{t, + 3p) (61) 

-'^--¥-^1 (62) 

= -3(«2 + /32)(/i + 3p + 6a2/32) (63) 

dy 

v^=a2{^i+p). (64) 



dy 

By introducing B = b^, = these equations can also be written as 

= -25/32 (65) 

dy 

^d(h ^ 5 _ 3^ (gg) 

«y 

^;^ = 02(02 +/32)-f^ (67) 
dy 

= ^2n{a2 + 2P2) (68) 

dy 

The equation of state fi — fi{p) is given by a lengthy third-order nonlinear differential 
equation, /i = constant is the interior Schwarzschild solution (see Theorem 1). 

We see that in the static limit uj2 0, we obtain Hi = H2 = and ei/i = e2/x — 0. Thus 
in this limit we obtain the interior Schwarzschild solution. The acceleration 02 is expressible 
in terms of p: 

al^j^i + p + DOi+p)\ (69) 

where D is a constant. For spherical symmetry, the quantities in the set S are invariant under 
the action of a 3-dimensional isotropy group at the center of the spheres. Thus, denoting 
center values by an index c, (a2)c = 0,^ and 



(70) 

3 (/ic +Pc^^ 



We collect our results in the following theorem: 

Theorem 2. The only rigidly rotating axistationary perfect fluid solutions with a pure mag- 
netic Weyl tensor are the interior Schwarzschild solution and the space-time given by Eqs. 



4 Restoring the metric 

We can always generate the set S from the metric (|l|) in the boosted frame (^. Thus, to 
obtain the full metric we only need to apply the commutators to the two coordinates t and 

•^This follows since a2 is the 2-component of the acceleration covector which must be zero at the center, 
otherwise there is a preferred direction. 
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of which all geometric quantities are independent, since the 1-forms lj^ and u>^ are known 
from solving the necessary and sufficient equations. Thus, to obtain the "missing part" of the 
metric, we must integrate 

[e,, Bj] t = {-t% ~ 7%) ent + {j^j, - j%) est (71) 
[e,, Bj] If = (7",-, - 7%) e„ip + (j^ji - 7%-) e^ip , (72) 

where the form of the right hand side follows from the fact that ei and 62 are linear combi- 



nations of dx and dy. Because of the complicated form of e^, the commutator equations (71) 



and (72) will in general be difficult to integrate. This is of no physical importance, though, 
since we know that the set S gives a complete local description of the space-time. The central 
dilemma is the condition for matching a given fluid solution to a vacuum counterpart. This 
problem can in principle be solved using the metric description, while no such procedure exists 
(yet) for the description in terms of S. 

4-.1 The shear- free case 

When the fluid is shear-free, i.e., rigidly rotating, we may put the boost parameter x equal 
to zero, since the frame (||) is comoving in the rigid case. The expressions ( p^ ) and ( [Til ) for 
Co and 63, respectively, simplify considerably: 

eo = r'^^dt (73) 

63 = — {d^ -todt) . (74) 
Q 

Here /, g and w are functions of x and y. 

We proceed to determining the x (angular) dependence of the metric quantities by use 
of the commutator equations (^ and (|7^). From the (0, 1) commutator acting on t, we get 
df /dx = 0, and hence / is a function of y only. Integrating the equation obtained from the 
(1,3) commutator acting on ip, we get 

g = giocos(^x) , (75) 



while acting on t gives 



-7^'^2£'osin(^x) +ujQ , (76) 
bf C 



where go and ujq depend on y only. The commutator equation (0,2) acting on t and (1,2) 
acting on x give, respectively 

(3^ = -y±lnC. (78) 
dy 

Comparison with (|6^) shows b = C( where one can set the constant C = 1 by rescaling the 
coordinate x. The action of (2, 3) on (p yields 

C = C— (79) 
Po 
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where we can arrange the constant C' — 1 by rescahng the coordinate ip. 

Substituting /?2 in the (2, 3) commutator acting on t , the resulting equation can be 
separated into x-dependent and a:-independent terms. From the x-independent part it follows 
that Wo is a constant. Hence we can set ujq — hy t ^ t + uiQip. The terms proportional to 
sin(5x) can be integrated, 

2uj2 - i^/fC' (80) 

where if is a constant. Putting together these results, we find uj — Ks'mx. By use of the 
transformation (p5|), we arrange v = \/J. The space-time metric takes the form (c/. Eq. (^) 

ds"^ = f{dt + Ksinxdip)^ ~ j[dy^ + plidx^ +cos^xd(p^)] . (81) 

The angular dependence of this metric is identical with the angular dependence of the 
metric in Ref. js). In that reference, the Ansatz is used that the eigenrays of the Killing 
vector d/dt are geodesic. Both solutions are of NUT type. For a full understanding of the 
relationship of the physical assumptions used in both papers, a study of the radial behavior 
is needed. 



5 Petrov types 

In this section, we classify the circularly rotating perfect fluid space-times into Petrov types. 
A Newman-Penrose null tetrad convenient for this investigation is given by 

l=l={eo + e3) (82) 

n =-^(60-63) (83) 

m = ^(ei +ie2) . (84) 
v2 

The nonvanishing Weyl spinor components are 

^-0 = ^ [E2 -E,- 2i/3 + iiH,-H2- 2E3)] (85) 
^2 = l[Ei+E2-iiH,+H2)] (86) 
*4 = ^ [E2 -E^ + 2i/3 + i{H^-H2 + 2E^)] . (87) 

The algebraic Petrov type is determined by the roots of the equation 

*4^''' + 6*2^^ + *o = . (88) 

The Petrov type is II if and only if either = or ^"4 = 0. (Reversal of the 63 direction 
interchanges and ^'4.) Petrov D type fields are characterized by the conditions 

AEiHi + 5E1H2 + 5E2H1 + AE2H2 - 2£'3iJ3 = (89) 
2El + bEiE2 + 2El -El- 2H^ - SiJjiJa - 2iJ| + iJ| = . (90) 

The type is N if ^'o = ^'2 = 0. Type III is excluded by equation (p8|). 
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From the point of view of star modehng, type N perfect fluids can hardly be of any sig- 
nificance. This is due to the following 

Theorem 3. A type N axistationary perfect fluid spacetime has a constant pressure p and 
density fi. 

Proof. The equations = ^2 ^ 0, with (||) and (||) yield E2 = H3 = -Ei and 
E3 — —H2 = Hi. We next arrange uj2 + (^2 — hy using the transformation (p^. We 
obtain Hi and Ei from the algebraic equations (^5|) and (p7|), respectively. The equation 
of state fi + p = follows at once from (|3C|). From the Bianchi identities ( ^6| ) and ( ^7| ) the 
statement of the theorem follows. □ 



5.1 Rigidly rotating electric space-times 

In this subsection we consider rigidly rotating fluid spacetimes with purely electric Weyl 
tensor. We then have cti = a2 = and Ha = 0. These spacetimes have the property stated 
in the following 

Theorem 4. A rigidly rotating purely electric perfect fluid space-time is of Petrov type D. 
The four-velocity lies in the plane spanned by the principal null directions of the Weyl tensor. 

Proof. We set cji = by using the transformation (|l^). From the Bianchi equations (p2[), 
and p5| ) we get, respectively 

E3 = (91) 

E2 = -2Ei (92) 

H = 'SE2+p. (93) 

The combination of Bianchi equations (^)-(pl|) yields 

a2 = (94) 
/?2 = -ai . (95) 



From the Ricci equations (pa) and (p4) we have, respectively 



ai = a2 (96) 
Ei^^iai(3i-al-iul) + ^p . (97) 



From conditions (89) and ( |90| ) it follows immediately that these fields are of type D. The 
double principal null directions are 1 + n + i(m — m) and 1 + n — i(m — m). Hence the four- 
velocity Co lies in the plane of the principal null directions. □ 

The condition that the rigidly rotating perfect fluid space-time is Petrov type D and that the 
four-velocity is a linear combination of the principal directions has been thoroughly investi- 
gated by Senovilla|9|, who gave all the exact solutions in this class. 
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